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t4.  Abstract  ^ 

;  ■  This  paper  reviews  radioactive  tracer  data  to  determine  effective  stratosphere- 
to- troposphere  transfer  rates  for  a  one-dimensional  parameterization.  The  basic  data 
used  are  the  time-dependent  global  stratospheric  burdens  of  95Zj) following  six  Chinese^ 
thermonuclear  explosions  Between^ 1 96 1976,  and  also  the  falloff  in  burdens  of C  Sr 
and  excess^v  C  following  the  very  large~Soviet  and  U.S.  test  series  of  1961  <aad^962. 

-Thp>jesults  of  the  Chinese  explosions  demonstrate  that  transfer  of  material  from 
stratosphere  to  troposphere  is  fast  in  winter  but  slow  in  the  summer  of  initial  in¬ 
jection.  It  is  best  modeled  by  a  seasonally  varying  eddy  diffusivity  profile,  show¬ 
ing  ‘'fast*  transport  (e.g.,  Dickinson-Chang,  Danielsen)  for  a  winter  injection,  and 
fast  transport  delayed  to  the  next  winter  for  a  summer  injection. 


The  cloud  rise  height  for  the  1961-1962  tests  is  not  well  known,  and  so  two  models 
for  the  yfeld-altitude  curve  are  used ; -following  Chang  et  al.  (1979)-;  > a  “high*  (Foley- 
Ruderman)  and  a  '"low1*  (Seitz)  model.  The  overall  time  variation^in  stratospheric  bur¬ 
den  of  and  in  stratospheric  and  tropospheric  burdens  of  can  be  simulated  best 

by  relatively  fast  falloff,  i.e.,  by  using  a  low  injection  height  (Seitz)  and  fast 
transport  (e.g.,  Dickinson-Chang).  The  difference  in  falloff  rate  of  stratospheric 
burdens  of  90g£\nd  0f  ^C  .is  due  principally  to  the  recycling  of  l^C^as  CO-J)  from 
rhe^troposphere ,  with  minor  corrections  due  to  sedimentation  of  the  aerosols  carrying 
90Srand  the  difference  in  effective  injection  height  of  these  two  tracers®  ‘ 90St  is 
a  tracer  of  fission  yield,  so  that  it  weights  the  smaller  bombs  with  a  low  ihjection 
height  more  heavily  than  does  ^C,  which  is  proportional  to  the  total  yield. 
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ABSTRACT 

This  paper  reviews  radioactive  tracer  data  to  determine 
effective  stratosphere-to-troposphere  transfer  rates  for  a 
one-dimensional  parameterization.  The  basic  data  used  are  the 
time-dependent  global  stratospheric  burdens  of  ^Zr  following 
six  Chinese  thermonuclear  explosions  between  1967  and  1976, 

90  ill 

and  also  the  falloff  in  burdens  of  Sr  and  excess  C  follow¬ 
ing  the  very  large  Soviet  and  U.S.  test  series  of  1961  and  1962. 

The  results  of  the  Chinese  explosions  demonstrate  that 
transfer  of  material  from  stratosphere  to  troposphere  is  fast 
in  winter  but  slow  in  the  summer  of  initial  injection.  It  is 
best  modeled  by  a  seasonally  varying  eddy  diffusivity  profile, 
showing  "fast”  transport  (e.g.,  Dickinson-Chang,  Danielsen) 
for  a  winter  injection,  and  fast  transport  delayed  to  the  next 
winter  for  a  summer  injection. 

The  cloud  rise  height  for  the  1961-1962  tests  is  not  well 

known,  and  so  two  models  for  the  yield-altitude  curve  are  used, 

following  Chang  et  al.  (1979):  a  "high"  (Poley-Ruderman)  and 

a  "low"  (Seitz)  model.  The  overall  time  variation  in  strato- 

90 

spheric  burden  of  Sr  and  in  stratospheric  and  tropospheric 
14 

burdens  of  C  can  be  simulated  best  by  relatively  fast  falloff, 
i.e.,  by  using  a  low  injection  height  (Seitz)  and  fast  trans¬ 
port  (e.g.,  Dickinson-Chang).  The  difference  in  falloff  rate 

90  14 

of.  stratospheric  burdens  of  Sr  and  of  C  is  due  principally 

lli 

to  the  recycling  of  C  (as  C02)  from  the  troposphere,  with 

minor  corrections  due  to  sedimentation  of  the  aerosols  carrying 
90 

Sr  and  the  difference  in  effective  injection  height  of  these 
90 

two  tracers:  Sr  is  a  tracer  of  fission  yield,  so  that  it 

weights  the  smaller  bombs  with  a  low  injection  height  more 

14 

heavily  than  does  C,  which  is  proportional  to  the  total  yield. 
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SUMMARY 


This  paper  reviews  radioactive  tracer  data  to  determine 
effective  stratosphere-to-troposphere  transfer  rates  in  a  one¬ 
dimensional  (1-D)  parameterization.  Atmospheric  motions  can¬ 
not  be  described  adequately  in  a  1-D  parameterization,  but 
1-D  codes  have  been  used  extensively  in  estimating  the  environ¬ 
mental  impact  of  atmospheric  lnjectants  such  as  oxides  of  nitro¬ 
gen  from  SSTs  and  halocarbons.  Thus  it  is  appropriate  to  eval¬ 
uate  different  1-D  transport  coefficients  (K-profiles ) ,  in  par¬ 
ticular  their  limitations. 

The  basic  data  used  here  are: 

•  The  time-dependent  global  stratospheric  burden  of  ^^Zr 
following  six  Chinese  thermonuclear  explosions  detonated 
at  Lop  Nor  (40°)  between  1967  and  1976. 

•  The  time-dependent  global  stratospheric  burdens  of  ^Sr 

14 

and  excess  C  following  the  very  large  Soviet  (polar) 
and  U.S.  (tropical)  nuclear  test  series  of  1961  and  1962. 

•  Ground-based  measurements  of  HTO  (T  =  %)  in  the  1960s. 
There  are  also  airborne  observations  of  HTO  from  1975 
onward,  but  the  data  are  still  being  reduced  and  thus 
are  not  included  here. 

Some  discussion  is  also  given  of  the  change  in  altitude  profiles 
of  9°Sr  and  li*C  following  the  Soviet  and  U.S.  test  series. 

These  tracers  all  have  somewhat  different  characteristics 
(Table  1-1,  p.  1-3)*  ^Zr  and  ^^Sr  are  fission  fragments,  so 
their  Intensity  is  proportional  to  fission  rather  than  total 
yield.  They  are  lodged  on  aerosols  which  are  subject  to  set¬ 
tling  and  are  precipitation  scavenged  in  the  troposphere. 

S-l 


HTO  is  a  gaseous  tracer  of  fusion  yield  which  is  precipitation 

scavenged  and  behaves  somewhat  like  water.  Excess  C,  which 

is  produced  by  reaction  of  neutrons  with  atmospheric  nitrogen, 

1^ 

is  a  tracer  of  total  yield.  It  is  carried  as  C02,  which  is 
a  good  tracer  for  atmospheric  C02.  Thus  its  fate  is  of  par¬ 
ticular  interest  from  a  climatic  standpoint.  CO,,  is  a  gas 
which  is  only  slightly  scavenged  by  precipitation  (unlike  the 
other  tracers)  and  is  only  lost  slowly  at  the  ground,  but  in 
contrast  to  the  other  tracers  it  is  recycled  from  the  tropo¬ 
sphere  back  into  the  stratosphere. 

First,  the  data  from  six  Chinese  thermonuclear  explosions 

are  analyzed.  These  explosions  had  very  similar  yields  of 

2-4  megatons  (Mt )  and  thus,  presumably,  similar  injection 

heights  ( —  18  km),  all  at  the  same  place  (Lop  Nor,  40°N, 

90°E),  which  provides  a  reasonable  if  not  ideal  simulation  for 

mid-latitude  injections  such  as  those  due  to  stratospheric 

aircraft,  whose  flight  path  peaks  near  50-55°N.  There  are  good 
95 

Zr  burden  data  from  all  six  explosions.  These  data  show 
rough  consistency  in  their  seasonal  variation,  with  slow  fall¬ 
out  in  summer  and  rapid  fallout  in  winter  (Fig.  2-2,  p.  2-6). 
There  are  also  limited  stratospheric  HTO  data  from  two  of  the 
six  Chinese  explosions.  These  data  are  still  being  reduced  and 
are  thus  excluded  here;  a  preliminary  examination  suggests  that 
the  HTO  data  are  roughly  consistent  with  the  ^Zr  data. 

For  winter  injections  a  "fast"  K-profile  (i.e.,  one  giving 
rapid  transport  near  the  tropopause,  such  as  DC  (Dickinson- 
Chang)  or  KD  (Danielsen)  (Fig.  4-1,  p.  4-2)  gives  better  re¬ 
sults  than  a  "slow"  K-profile,  such  as  WT  (Whitten  and  Turco 
or  "Wofsy-Type" ) ,  HN  (Hunten),  or  CH-LO  (Crutzen-Howard/Low) . 
For  times  of  one  to  three  years  after  a  summer  or  fall  injec¬ 
tion  it  is  best  to  use  a  "fast"  profile  with  the  injection 
delayed  to  the  following  winter,  or  possibly  a  "slow"  K-pro¬ 
file,  but  this  latter  does  not  appear  to  be  satisfactory  during 
the  initial  summer/fall  season  because  it  predicts  too  much 
initial  fallout. 


S-2 


In  view  of  the  potential  importance  of  the  seasonal  vari¬ 
ation  in  the  removal  of  injectants  from  the  stratosphere,  this 

95 

is  discussed  next,  and  both  the  Chinese  Zr  and  the  earlier 
U.S.  and  USSR  ^°Sr  fallout  are  reviewed.  For  summer,  fall, 
and  winter  injections  (there  are  no  data  on  springtime  injec¬ 
tions)  there  is  evidence  for  a  seasonal  variation,  with  very 
little  fallout  until  winter  (i.e.,  up  to  6  months  following  a 
summer  injection),  and  then  there  is  roughly  a  4-month  delay 

between  the  beginning  of  fallout  from  the  stratosphere  (January- 

90 

February)  and  the  peak  in  Sr  in  rain  at  the  ground,  which 
occurs  in  May-June. 

90 

Next,  the  time  variations  in  stratospheric  burdens  of  7  Sr 

and  (as  ^COg)  following  the  1961-1962  test  series  are 

analyzed.  Figure  2-7  (p.  2-13)  shows  the  changing  stratospheric 

burdens  of  these  tracers  (and  of  HTO)  and  the  corresponding 

14 

tropospheric  burden  of  excess  CO^.  The  stratospheric  burden 
of  9°Sr  decreases  much  faster  than  that  of  1^C,  with  the  differ¬ 
ence  maximizing  some  2-4  years  after  injection.  The  tropospheric 

burden  of  li+C,  some  of  which  comes  from  the  injection  of  the 

90 

1950s,  is  large;  by  contrast,  the  tropospheric  burden  of  '  Sr 
is  very  small. 

Before  one  can  model  the  time  dependence  of  atmospheric 
burdens  of  and  ^°Sr,  one  must  determine  the  effective  in¬ 
jection  height,  because  in  the  1950s  and  1960s  there  were  a 
variety  of  injections  at  different  latitudes  (U.S.  injections 
equatorial,  Soviet  injections  mainly  polar)  and  altitudes  or 
yields.  Two  different  models  of  cloud  height  as  a  function  of 
yield,  as  used  by  Chang  et  al.  (1979) >  are  considered  here. 

A  review  of  data  finds  little  evidence  for  the  higher  cloud- 
rise  model. 


S-3 


on 

Computer  modeling  of  the  change  in  stratospheric  Sr  and 
14 

excess  C  burdens,  using  all  the  variations  discussed  here, 
gives  the  following  results: 

1.  ^Sr  falls  out  faster  than  ^C,  principally  because  of 
its  different  behavior  near  and  on  the  ground,  al¬ 
though  there  are  small  contributions  due  to  the  sedi- 
mentation  of  Sr  and  its  lower  effective  injection 
height  (it  is  a  tracer  of  fission  yield,  while  C  is 
a  tracer  of  total  yield:  the  biggest  bombs,  which  go 
to  the  greatest  heights,  have  a  relatively  small  fis¬ 
sion  fraction,  so  that  ^Sr  in  1961-1962  was  injected 
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1-2  km  lower,  on  the  average,  than  C). 

2.  The  decline  in  burdens  (stratospheric  and  tropospheric 

for  and  stratospheric  for  ^°Sr)  is  best  explained 

by  assuming  a  low  injection  height  (which  is  in  accord 
with  most  available  data)  and  fast  transport  (which 
agrees  with  the  conclusion  for  the  Chinese  tests). 

See  Pig.  6-2  (p .  6-4)  for  ^°Sr  and  Pig.  6-3  (pp.  6-7, 
6-8)  for  14C02. 

3.  There  is  no  need  to  postulate  excessively  rapid  fall¬ 
out  of  9°Sr,  or  a  large  reservoir  for  ^C  at  high 
altitudes . 

Johnston,  Kattenhorn,  and  Whitten  (1976)  studied  the 
change  in  altitude  profile  of  ^C  between  1963  and  1964-1970 
using  a  1-D  representation  of  the  data,  and  found  that  a  rela¬ 
tively  slow  transport  profile  (HN,  i.e.,  Hunten)  best  fitted 
the  data.  The  work  of  Johnston  et  al.  has  been  reviewed  care¬ 
fully,  since  their  conclusion  differs  from  that  of  the  present 
study,  and  the  following  conclusions  are  drawn: 

1.  The  analysis  of  Johnston  et  al.  uses  boundary  con¬ 
ditions  which  we  believe  to  be  inappropriate  for  ^C, 
with  a  strong  sink  at  the  bottom  (1  km)  and  another 
sink  at  the  top  (50  km). 


2.  Since  no  1-D  parameterization  can  be  expected  to  match 
the  full  complexity  of  stratosphere-to-troposphere  trans¬ 
fer,  should  one  intercompare  the  altitude  profiles  of 
mixing  ratio  or  of  number  density  of  tracer,  or  the 
total  stratospheric/tropospheric  burdens,  and  at  what 
times  after  injection?  Johnston  et  al.  emphasize 
matching  profiles  of  mixing  ratio  and  of  number  density 
at  relatively  short  intervals  by  re-initializing  their 
computation  (so  that  the  effect  of  the  boundary  con¬ 
ditions  and  of  interhemispheric  transfer  is  minimized). 

By  contrast,  the  present  work  emphasizes  the  total 
stratospheric  burden  at  relatively  long  times  after 

inj  ection. 

3.  Explicitly,  2-4  years  after  an  effective  injection  in 

January  1963,  one  can  get  the  same  stratospheric  burden 
14 

of  C  either  by  using  fast  (DC)  transport  with  what 
we  believe  to  be  the  correct  lower  boundary  condition 
or  by  using  relatively  slow  (HN)  transport  with  an  ex¬ 
cessively  strong  sink  at  the  bottom  (Fig.. 6-5,  p.  6-13). 

In  conclusion,  the  present  analysis  provides  a  consistent 

explanation  of  available  data  for  the  decay  in  stratospheric 

95  90 

burdens  of  Zr  following  Chinese  explosions,  and  of  Sr  and 

14 

C  following  U.S.  and  USSR  tests.  Of  the  various  factors  con¬ 
sidered  here,  the  following  are  found  to  be  in  agreement  with 
the  data: 

1.  Seasonal  variation  in  initial  fallout  from  Chinese 
tests,  with  little  or  no  fallout  following  a  summer  or 
fall  injection  until  the  following  winter  season 

2.  Fast  transport  (i.e.,  DC  or  KD  model,  rather  than  WT) 
following  the  1961-1962  tests,  and  following  Chinese 
injections  In  winter 

3.  Low  injection-height  models 
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4.  Generally  accepted  physics  for  the  boundary  conditions 

of  different  tracers,  for  the  sedimentation  of  aerosols, 

qq  1 4 

and  for  the  relative  injection  height  of  Sr  and  C. 

However,  no  1-D  eddy  diffuslvlty  profile  can  describe  the  com¬ 
plexity  of  actual  transport  over  all  relevant  time  scales,  and 
so  the  user  is  urged  to  be  sure  that  the  parameterization  is 
appropriate  to  the  problem  at  hand. 


1.  INTRODUCTION 


In  investigating  the  stratospheric  pollution  due  to  air¬ 
craft  and  other  artificial  sources,  one  key  problem  has  to  do 
with  how  long  injectants  remain  in  the  stratosphere,  since  for 
chemically  active  injectants  such  as  N0X  (=  N0+N02+...), 

H0x  (=  OH+HC^+I^Op ) ,  and  C1X  (=  C1+C10+...)  this  is  one  crucial 
factor  in  determining  their  effect  on  the  overall  stratospheric 
photochemistry.  At  present,  much  or  most  of  the  modeling  used 
for  making  policy  decisions  on  stratospheric  pollution  treats 
the  motions  one-dimensionally ,  because  the  one-dimensional 
(1-D)  models  can  be  run  more  quickly  and  used  more  economically 
In  sensitivity  studies  than  can  more  complex,  higher  dimensional 
models . 

It  is  evidently  important  to  validate  the  transport  models 
used.  They  have  generally  been  based  on  the  observed  profiles 
of  molecules  such  as  N20  and  CH^ ,  which  have  a  steady  ground- 
level  source,  a  fairly  long  atmospheric  lifetime,  and  a  postu¬ 
lated  photochemical  loss  mechanism  In  the  upper  stratosphere. 

By  contrast,  radioactive  tracers  due  to  thermonuclear  explosions 
provide  an  instantaneous  point  injection,  generally  present  no 
chemistry,  and  have  a  relatively  simple  loss  mechanism  by  wet 
or  dry  deposition  on  the  ground  or  in  the  ocean.  (Note  that 

llj  il} 

radioactive  C  Is  transformed  rapidly  into  C02  and  then 
follows  this  species.) 

The  principal  source  of  data  for  all  radioactive  tracers 
comes  from  the  U.S.  Department  of  Energy  "Airstream"  flights 
and  the  preceding  Defense  Nuclear  Agency  "Stardust"  flights, 
which  for  the  past  25  years  have  operated  several  times  per 
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year  over  a  corridor  near  the  west  coast  of  America  from  75°N 
to  the  equator  (to  50°S  through  April  1974)  at  altitudes  up  to 
19.2  km.  These  aircraft  data  have  been  supplemented  by  balloon 
flights,  mainly  near  30°N  and  30°S. 

Here  I  first  examine  the  falloff  with  time  of  the  strato- 
95 

spheric  burden  of  Zr  due  to  six  Chinese  thermonuclear  explo¬ 
sions  between  June  1967  and  November  1976.  The  same  technique 

90  111 

is  then  applied  to  the  stratospheric  burdens  of  Sr  and  C  due 
to  the  Soviet  and  U.S.  test  series  of  1961-1962.  The  charac¬ 
teristics  of  the  different  tracers  are  given  in  Table  1-1,  and 
those  of  the  Chinese  versus  Soviet  and  U.S.  series  are  given 
in  Table  1-2.  The  Soviet  series  provided  a  very  large  source 
(in  both  space  and  time),  at  a  period  when  there  were  extensive 
measurement  networks  operational  on  a  worldwide  basis.  By  con¬ 
trast,  the  Chinese  tests  provided  fairly  reproducible  point  in¬ 
jections  of  2-4  Mt  yield  at  a  time  when  only  limited  fallout 
measurements  were  made.  Regarding  the  tracers,  Sr  and  ^Zr 

are  physically  comparable  tracers,  precipitation-scavengeable 

95 

fission  fragments.  Zr  has  a  short  half-life,  i.e.,  a  high 
initial  activity  per  unit  mass  but  correspondingly  a  short 
available  observation  time  ( ^  1  year),  while  Sr  has  a  longer 
half-life,  a  lower  initial  activity  per  unit  mass,  and  a  cor¬ 
respondingly  longer  available  observation  time.  The  concen¬ 
tration  of  14C  is  proportional  to  total  yield  (Mt)  of  the  bomb, 

that  of  T  is  proportional  to  fusion  yield,  and  those  of  ^®Sr  and 
95  lit 

^Zr  are  proportional  to  fission  yield.  C  is  only  slightly 

scavenged  by  precipitation;  it  is  a  gaseous  tracer  vrhich  is  re- 

14 

cycled  from  the  troposphere  (as  C02)  and  thus  provides  quite 
different  information  than  do  the  fission  fragments  or  HTO. 

HTO  is  scavenged  by  precipitation  and  is  also  lost  by  vapor 
exchange,  so  that  it  is  considered  to  be  a  tracer  for  water. 

For  an  application  to  stratospheric  aircraft  flight,  all 
the  tracers  have  the  advantage  that  they  are  stratospheric 
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TABLE  1-1.  CHARACTERISTICS  OF  DIFFERENT  TRACERS 
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Some  atmospheric  HTO  data  from  1975 


point  injections  in  approximately  the  correct  altitude  range 
(for  SSTs),  around  18  to  20  km;  the  Chinese  mid-latitude  injec¬ 
tions  provide  the  closest  simulation  of  the  latitudes  of  air¬ 
craft  (approximately  50-55°N)  and  also  provide  seasonal  varia¬ 
bility,  while  the  Soviet  high-latitude  injections  provide  a 
higher-altitude,  high-latitude  source  and  a  measure  of  recyc¬ 
ling  from  the  lower  atmosphere,  i.e.,  longer  time  behavior  than 
can  be  observed  from  the  lower-yield  Chinese  explosions. 

Table  1-3  amplifies  this  discussion  by  indicating  what 
data  are  used  and  what  information  is  gathered  from  a  given 
set  of  data  so  as  to  be  able  to  decide  between  different  K- 
profiles,  which  represent  the  customary  1-D  parameterization 
of  dynamics.  Section  2  reviews  the  data  used.  Section  3  treats 
the  problem  of  cloud  rise  height.  Section  ^  discusses  different 
K-profiles,  and  Section  5  indicates  details  of  the  1-D  diffusion 
model  used  here.  Comparison  between  data  and  computations  is 
made  in  Section  6,  and  conclusions  are  presented  in  Section  7. 

It  must  be  stressed  that  the  atmosphere  Is  not  one-dimen¬ 
sional,  and  thus  any  1-D  parameterization  Is  very  limited.  In 
particular: 

•  1-D  models  and  correspondingly  K-profiles  are  Inter¬ 
preted  In  three  distinct  ways  that  are  not  always 
compatible : 

1.  As  representing  atmospheric  conditions  near  30°N, 
or  generally  at  mid-latitude 

2.  As  representing  a  hemispheric  average 

3.  As  representing  a  global  average  for  a  Northern 
Hemisphere  mid-latitude  injection. 

•  It  will  be  shown  that  no  single  K-profile  can  account 
for  all  aspects  of  available  radioactive  tracer  data 
(e.g.,  burden  and  profile  at  all  times  after  injection). 
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but  nevertheless  these  limitations  show  deficiencies 
of  the  1-D  parameterization  rather  than  evidence  of 
a  lack  of  physical  understanding. 


TABLE  1-3.  DATA  USED  AND  INFORMATION  GATHERED 


a.  Data  Used 

a.l.  Zr  1967-78  stratospheric  burden  from  Telegadas  1974,  1976,  1979b. 

a. 2.  HTO  1963-71  ground  level  data,  IAEA  (see  especially  Krey  and 
Krajewski  1970). 

90 

a. 3.  Sr  1954-72  stratospheric  burden;  data  from  National  Academy  of 
Sciences  1975,  p.  110;  data  verified  from  Friend  et  al.  1961, 

Feely  and  Troutman  1971,  and  Leifer  and  Toonkel  1978. 

14 

a. 4.  C  1955-69  stratospheric  burden;  data  from  Telegadas  1971. 

a.  5.  90Sr  and  14C  altitude  profile  of  Johnston  et  al.  1976. 

b.  Information  Sought 

b. l.  Overall  stratosphere-to-troposphere  transport,  including  seasonal 

variability,  from  95zr  data  (a.l). 

90 

b.2.  Testing  of  K-profiles,  from  Sr  data  (a. 3). 

14  90 

b. 3.  Short-term  variability  between  C,  HTO,  and  Sr  may  tell  us  about 

sedimentation,  effective  injection  height  differences,  and  differ¬ 
ences  in  boundary  conditions. 

c.  Ultimate  Results 

c. l.  Intercomparison  of  different  K-profiles  in  stratosphere  and  perhaps 

also  in  troposphere. 

c.2.  Verification  of  the  boundary  conditions. 

c.3.  Some  sketchy  information  on  injection  height  of  bombs. 
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2.  DATA  USED  HERE 


2.1 .  95Zr  DATA 

Telegadas  (1974,  1976,  1979c, d)  provides  altitude-latitude 
profiles  and  also  total  stratospheric  burdens  of  ^Zr  from  1967 
to  1978.  Pig.  2-1  shows  the  stratospheric  burden  of  9^Zr 
(decay-corrected  to  the  date  of  release)  for  the  six  Chinese 
thermonuclear  explosions  that  have  been  detonated  in  the  at¬ 
mosphere.  By  far  the  largest  part  of  the  burden  lies  in  the 
Northern  Hemisphere  (see,  e.g.,  Telegadas  1974,  Pigs.  47-49). 

Table  2-1  lists  some  characteristics  of  the  Chinese  at¬ 
mospheric  thermonuclear  explosions  considered  here.  Note  that 

there  were  also  some  French  thermonuclear  explosions  which 

95 

contributed  to  the  Zr  burden,  but  these  have  been  separated 
from  the  Chinese  debris  in  Telegadas'  analysis  of  the  data. 
Table  2-2  lists  the  actual  data  that  are  used  here. 

From  Fig.  2-1  one  sees  that  the  stratospheric  burden  re¬ 
mains  constant  after  a  summer  or  fall  injection  until  the 
following  January  or  February,  when  significant  fallout  of  the 
debris  begins.  The  actual  date  varies  a  little  from  year  to 
year;  Telegadas'  estimates  are  listed  in  Table  2-3.  To  display 
this  seasonal  variability  of  the  data,  in  Fig.  2-2  are  shown 
the  stratospheric  burdens  normalized  [to  the  maximum  burden 
injected  (Table  2-1)]  and  plotted  on  a  seasonally  adjusted 
basis. 

In  this  way  all  the  data  fall  more  or  less  on  top  of  one 

95 

another,  showing  a  synthetic  1-  to  2-year  model  of  how  Zr, 
or  presumably  any  radioactive  aerosol,  falls  out  of  the  strato¬ 
sphere  after  a  pulsed  injection  at  mid-latitude  (40°N)  and 
approximately  18  km  altitude  (Telegadas  1974,  Figs.  44-45). 
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FIGURE  2-1.  g5 

Stratospheric  yDZr  burden  (decay 
corrected)  following  the  Chinese 
nuclear  tests.  Observed  burdens 
are  indicated  by  the  short  hori¬ 
zontal  lines.  Note  that,  regard¬ 
less  of  the  date  of  injection, 
significant  fallout  from  the 
stratosphere  starts  in  the  North¬ 
ern  Hemisphere  winter  season,  de¬ 
monstrating  that  stratosphere-to- 
troposphere  transfer  is  much 
faster  in  winter  than  in  summer 
of  the  injection  hemisphere. 
(Source:  Telegadas  1974,  1976, 
1979b) 
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TABLE  2-1.  LIST  OF  CHINESE  ATMOSPHERIC  THERMONUCLEAR 
EXPLOSIONS  (ALL  AT  LOP  NOR,  40ON,  9QOE) 


TABLE  2-2.  STRATOSPHERIC  INVENTORY  OF  ZIRCONIUM-95 
RESULTING  FROM  CHINESE  TESTS 

(Source:  Telecjadas  1  974  ,  1  976,  1  979c,  d } 


Date 

Predominantly 

Chinese  Debris,  kCi 

17  June  1967  Test  ( decay-corrected  to  17  June  1967) 

October  7  -  November  9,  1967 

38,500 

January  8  -  January  21,  1968 

36,400 

April  1  -  April  7,  1968 

27,200 

May  30  -  June  15,  1968 

23,000 

27  December  1968  Test  (decay-corrected  to  27  December  1968) 

February  3  -  February  19,  1968 

38,500 

April  6  -  April  9,  1969 

27,800 

July  14  -  July  29,  1969 

19,500 

29  September  1969  Test  (decay-corrected  to  29  September  1969) 

October  13  -  October  17,  1969 

44,700 

January  5  -  January  20,  1970 

45,300 

May  6  -  May  11 ,  1970 

i  32,200 

14  October  1970  Test  (decay-corrected  to  14  October  1970) 

February  22  -  February  27,  1971 

38,400 

May  22  -  May  28,  1971 

29,800 

27  June  1973  Test  (decay-corrected  to  27  June  1973) 

September  5  -  September  19,  1973 

33,200 

October  30  -  November  9,  1973 

34,700 

January  22  -  February  2,  1974 

28,300 

April  12  -  April  29,  1974 

19,800 

17  November  1976  Test  (decay-corrected  to  17  November  1976) 

March  22  -  April  10,  1977 

58,100 

July  6-22,  1977 

45,500 

October  12-29,  1977 

36,100 

April  6-21  ,  1978 

24,700 

TABLE  2-3.  ESTIMATED  START  OF  SIGNIFICANT 
FALLOUT  FOLLOWING  CHINESE  EXPLOSIONS 

(Source:  K.  Telegadas,  personal  communication) 


Code 


M 

M 

M 

M 

M 

M 


1 

2 

3 

4 

5 

6 


Injection  Date 
17  Jun  67 
27  Dec  68 
29  Sep  69 
14  Oct  70 
27  Jun  73 
17  Nov  76 


Start  of  Fallout 
15  Jan  68 
15  Jan  69 
15  Feb  70 
15  Feb  71 
15  Dec  73 
1  Feb  77 
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For  comparison  with  Fig.  2-2,  which  demonstrates  the  sea¬ 
sonal  character  of  fallout,  Fig.  2-3  shows  the  same  normalized 
burdens  as  a  function  of  time  after  injection.  I  think  it  is 
reasonable  to  conclude,  as  did  Telegadas,  1974,  that  seasonal 
variation  is  a  significant  factor  in  the  fallout.  Further, 

Fig.  2-4  plots  the  delay  time  between  the  injection  (from 
Table  2-1)  and  the  start  of  significant  fallout  (from  Table 
2-3).  This  again  shows  a  definite  relation  between  the  month 
of  injection  and  the  time  T^  of  initial  fallout  that  is  con¬ 
sistent  with  the  observation  that  substantial  fallout  tends  to 
begin  in  the  spring  following  the  initial  injection. 

The  seasonal  variations  following  the  nuclear  tests  of 
the  1950s  and  1961-1962  are  discussed  later  (Section  2.5),  but 
first  I  discuss  HTO  data  following  the  Chinese  explosions. 

2.2.  HTO  DATA 

During  the  1960s  there  was  an  extensive  network  of  ground- 
based  stations  for  the  collection  of  radioactive  fallout.  One 

O 

tracer  of  interest  here  is  tritium  (JH  *  T),  which  is  a  product 
of  fusion.  It  is  carried  mainly  as  tritiated  water,  HTO,  which 
can  be  a  fairly  good  tracer  for  water  vapor  in  the  atmosphere. 

HTO  is  normally  considered  to  be  precipitation  scavenged 
in  the  troposphere,  and  thus  the  ground-based  fallout  is  taken 
as  a  measure  of  loss  from  the  stratosphere.  Following  the 
1961-1962  tests  there  was  evidence  of  a  decline  In  strato¬ 
spheric  burden  with  a  half-life  of  16  months  or  an  e-folding 
time  of  23  months  (Krey  &  Krajewski  1970;  see  also  IAEA  1969, 
etc . ) . 

Since  1975  airborne  measurements  of  HTO  and  HT  have  been 
made  a  part  of  the  Airstream  program  (Mason  1978).  It  had 
been  hoped  to  analyze  these  data  in  the  present  work,  but  this 
has  not  been  done  because  of  delays  in  the  resolution  of  some 
difficulties.  A  preliminary  examination  suggests  that  the  HTO 
data  are  generally  consistent  with  the  7  Zr  results  (personal 
communication,  A.  S.  Mason  &  K.  Telegadas). 
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FIGURE  2-4 .  Delay  time  T  i  between  date  of  injection  (Table  2-1) 
and  date  of  initial  fallout  (Table  2-3)  for  six 
Chinese  thermonuclear  explosions. 
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2.3.  STRATOSPHERIC  BURDENS  OF  90Sr  AND  14C  DURING  THE  1960s 

90 

Figure  2-5  shows  the  global  stratospheric  burden  of  Sr 
for  the  period  195^-1972  (National  Academy  of  Sciences  1975, 
p.  110,  Fig.  A. 19;  data  verified  from  Friend  et  al.  1961, 

Feely  &  Trautman  1971,  and  Leifer  &  Toonkel  1978),  and  Fig. 

.  i  ii 

2-6  shows  the  stratospheric  burden  of  C  for  the  period  1955 

to  1969  (National  Academy  of  Sciences,  op.  cit.,  verified  from 

Telegadas  1971).  For  a  comparative  analysis  the  stratospheric 
1*1  90 

burdens  of  C  and  Sr  are  normalized  to  the  same  (maximum) 

value  in  January  1963,  and  they  are  plotted  together  in  Fig. 

2-7,  which  also  shows  the  difference  between  the  (normalized) 

90  lA 

stratospheric  burdens  of  Sr  and  C  and  the  tropospheric 

i  ii 

burden  of  C  (from  Telegadas  1971).  Figure  2-7  also  shows 
the  stratospheric  burden  of  HTO  inferred  from  ground-based 
data  (Krey  &  Krajewski  1970). 

Figure  2-7  will  be  used  for  comparison  with  different 
numerical  models.  A  number  of  points  should  be  made  here: 

•  Table  2-4  shows  the  overall  magnitude  of  the  various 
atmospheric  thermonuclear  test  series:  the  1961-1962 
test  series  provided  the  biggest  input,  but  a  signifi¬ 
cant  background  remained  due  to  the  testing  of  the 
later  1950s.  One  sees  that  the  Chinese  tests  of  1967- 
1976  correspond  to  much  smaller  injections. 

•  Table  2-5  gives  a  detailed  breakdown  of  the  1961-1962 

test  series,  showing  both  fission  and  total  yields; 

90 

recall  that  Sr  production  is  proportional  to  fission 
14 

yield,  while  C  production  is  due  to  total  yield. 

•  On  the  scale  of  Fig.  2-7,  where  the  peak  injection 
corresponds  to  3^0  Mt  total  or  96  Mt  of  fission,  the 
peak  injection  due  to  a  typical  Chinese  bomb,  which 
has  1.5  Mt  fission  or  3  Mt  total  yield,  would  be 
3/340  —  0.01  for  14C,  1.5/244  a  0.006  for  HTO,  and 


2-10 


FIGURE  2-5.  Global  stratospheric  burden  of  ^Sr,  1954-1972. 
(Source:  Table  1-3) 
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TABLE  2-4.  CHRONOLOGY  OF  ATMOSPHERIC  NUCLEAR  TESTS 


Time 

Country/Location 

Total  Yield,  Mt 

1952-54 

U.S.  (Tropical ) 

59  l 

60 

USSR  (Mid-hiah  latitude) 

A 

1955-56 

U.S.  (Tropical) 

1  9  ( 

28 

i 

USSR  (Mid-high  latitude) 

9  j 

1957-58 

U.S.  (Tropical) 

30  j 

| 

U. K.  (Tropical ) 

7  , 

>  85 

USSR  (Mid-high  latitude) 

48  J 

1 

1959-60 

Morator i um : 

No  atmospheric  testing 

U 

1961 

USSR  (High  latitude) 

120 

120 

1962 

U.S.  (Tropical) 

37  i 

j  220 

USSR  (High  latitude) 

1 80  1 

1963 

Partial  test  ban  treaty: 

No  atmospheric  testing 

u 

1967-68 

France  (Tropical ) 

4  ( 

!  io 

i 

China  (Mid-latitude) 

6I 

1969-70 

France  (Tropical ) 

2  I 

! 

China  (Mid-latitude) 

6 ! 

1971 

France  (Tropical) 

i 

1 

1973 

China  (Mid-latitude) 

2-3 

2-3 

1976 

China  (Mid-latitude) 

4 

4 
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an  gc 

1.5/96  »  0.016  for  Sr  or  Zr,  which  represent,  re¬ 
spectively,  total  yield,  fusion  yield,  and  fission 
yield. 

•  Note  that  there  appears  to  be  a  problem  with  the  total 
14C  budget  (see  Pig.  2-8)  which  suggests  a  larger  in- 

14 

put  of  C  than  is  observed  with  the  global  measure¬ 
ments.  This  may  well  be  due  to  an  uncertainty  in  the 

14 

production  rate  or  the  number  of  C  atoms/Mt  (see 
Machta  et  al.  1963). 

•  I  shall  only  use  the  data  following  the  1961-1962 
test  series,  but  note  the  residual  effects  due  to 
the  earlier  tests,  which  must  be  taken  into  account 
as  providing  a  background. 

90  14 

•  While  the  stratospheric  burdens  of  Sr  and  C  show 
general  agreement,  there  are  a  number  of  differences: 

90 

-  The  stratospheric  Sr  burden  falls  off  much  more 

14 

rapidly  than  does  the  C  burden. 

14 

-  Note  the  large  C  burden  in  the  troposphere;  the 

90 

tropospheric  Sr  burden  is  very  small. 

•  Note  that  the  stratospheric  HTO  burden  falls  off  al- 

90 

most  as  rapidly  as  that  of  Sr  '’nd  definitely  faster 
14 

than  that  of  C,  especially  at  times  in  excess  of 
2.5  years. 

2.4.  SEASONAL  VARIATION  OF  FALLOUT  DURING  THE  1950s  AND  1960s 

To  supplement  the  discussion  of  seasonal  variation  in 
fallout  from  the  Chinese  tests  (see  Section  2.1,  in  particular 
Pig.  2-4),  let  us  ask  what  can  be  learned  from  ^°Sr  fallout 
data  following  the  major  thermonuclear  tests  of  the  1950s  and 
1960s.  A  good  data  base  for  this  is  provided  by  HASL  1977, 
which  is  a  tabulation  of  all  available  monthly  fallout  data 
of  ^Sr  from  1954  through  1976.  While  there  exist  monthly 
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1964 


1966  1968  1970 


1960  1962  1964  1956  1956  1960  1962 

U.S.  U.S.,  U  K.  U.S.,  U.K.  U.S.,  USSR  USSR  USSR 
USSR  USSR  U.K.  (U.S.) 

MM»ir 

FIGURE  2-8.  Excess  s t ra tospheri c  from  atmospheric  thermo¬ 
nuclear  explosions.  The  total  injected  is 
shown  [2  x  10‘°  atoms/Mt  for  airbursts,  half  that 
for  surface  bursts  (Machta  et  al .  1963)]  as  well 
as  the  origin  of  the  injection  (U.S.,  (J.K.,  F  = 
tropical,  USSR  =  arctic,  C  =  mi d- 1 ati tude ) . 
(Source  of  data:  Telegadas  et  al.  1971) 


global  data  from  1958  through  the  early  1970s  and  quarterly 

data  thereafter  (see,  e.g.,  Volchok  &  Kleinman  1971,  Feely 

1976,  Harley  1976).  I  shall  use  the  monthly  deposition  data 

from  New  York  City,  which  go  from  1954  to  the  present,  and 

which  provide  information  comparable  to  that  of  Fig.  2-4  for 

89  90 

the  Chinese  tests,  with  an  indication  (from  the  Sr/  Sr 
ratio)  of  how  old  the  debris  is. 

The  major  stratospheric  injections  are  listed  in  Table 

2-6,  with  an  indication  (from  HASL  1977,  p.  A-73ff . )  of  when 

(in  what  month)  the  peak  in  fallout  occurred.  The  annual 

peak  in  fallout  deposition  (maximum  monthly  rate  of  fallout, 

90  2 

MCi  Sr/km  )  occurs  generally  between  April  and  June,  whereas 
the  beginning  of  fallout  from  the  stratosphere  shown  in  Fig. 

2-2  generally  occurs  in  January-February .  Thus  there  is  a 
certain  time  delay,  typically  four  months,  between  the  be¬ 
ginning  and  the  maximum  of  fallout. 

In  Fig.  2-9  are  shown  the  time  delay  data  from  Table  2-6, 
including  the  data  for  the  Chinese  tests  from  Fig.  2-4,  plotted 
as  a  function  of  the  month  of  deposition.  Each  point  is  la¬ 
beled  with  a  letter  for  the  latitude  of  injection  [T  =  tropical, 
(i.e.,  U.S.),  M  =  mid-latitude  (i.e.,  China),  P  =  polar  (I.e., 
USSR)]  together  with  a  number  (see  Table  2-6)  to  characterize 
the  specific  Injection.  Points  P3  and  T4  correspond  to  injec¬ 
tions  over  rather  long  time  intervals  and  thus  should  not  be 
given  as  much  weight  as  the  others.  Nevertheless,  it  Is  clear 
that  the  data  plotted  In  this  way  show  a  significant  seasonal 
variation  In  fallout  deposition,  with  the  debris  from  injec¬ 
tions  at  all  latitudes  being  deposited  in  the  spring. 

Note  also  a  current  observation  by  Holloway,  1979 s  who 

finds  that  following  the  Chinese  June  1973  injection  the  (decay- 

v  91  144 

corrected)  surface  deposition  of  Y  and  Ce  in  July-October 

1973  was  1S5  of  that  in  March-June  1974. 
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TABLE  2-6.  DELAY  BETWEEN  STRATOSPHERIC  INJECTION  AND  FALLOUT 


A.  U.S.  4  Soviet  Tests.  1956-1962a 


Code 

Yield,  Mt/ 

Dates 

Max.  Surface 

Delay  T. 

NO  . 

Coun try/ Lat i tude  No.  Tests 

Mean 

Range 

Fallout  Rate 

months ‘ 

Ti 

U . S . /Trop i cal 

10.5/5 

6/56 

5-7 

4/57 

10tl 

?z 

USSR/Mi  d-high 

6/4 

9/56 

8-4 

4/57 

7+  1 
-2 

F  3 

USSR' Mid-high 

12/5 

6/  :  7 

4;  8-10 

5/58 

(11) 

T  4 

U.S.  4  U.K. /Tropical  21/30 

6/58 

11/57;  4-9 

4/59 

(10) 

P  5 

USSR/ Arcti c 

36/10 

10/58 

-- 

4/59 

6 

P6A 

USSR/Arcti c 

37/13 

10/61 

9-11 

2/62 

4 1 1 

p  6  B 

USSR/Arcti c 

83/2 

11/61 

-- 

7/63 

8 

T7A 

U . S . /T  ropi ca  1  j 

7/62 

4-7 

7/63 

12  +  3 

T7B 

U . S . / T ropi cal  \ 

37/7 

10/62 

10-11 

7/63 

6-1 

°S 

USSR/Arcti c 

180/16 

10/62 

8-12 

7/63 

9*2 

B. 

Chinese  Tests  (all 

at  40°N) 

Code 

Initial 

Fallout 

Delay  Times, 

months 

No. 

Date  Fal loutb 

Maximum^  to  Initial  Fallout  (T,)  to 

Fallout  Maximum  (T0) 

Ml 

6/67  15  Jan 

5/68 

7 

n 

M2 

12/68  15  Jan 

5/69 

0.5 

6 

M3 

9/69  15  Feb 

6/70 

4.5 

9 

M4 

10/70  15  Feb 

5/71 

4.0 

7 

M5 

6/73  15  Dec 

5/74 

5.5 

11 

M6 

11/76  1  Feb 

_  _  d 

2.5 

-- 

injections:  see,  e.g.,  Bauer  1978,  Table  B-7,  p.B-21ff. 
Fallout:  see  HASl  1977. 

bFrom  Telegadas,  personal  communication. 
cFrom  HASL  1977. 
d0nly  quarterly  data. 


DELAY  TIME  T2,  months 


Reference  should  be  made  to  earlier  discussions  of  NOAA 
(e.g.,  Machta  1965,  Telegadas  1974)  which  point  out  that  care¬ 
ful,  long-term  observers  of  fallout  conclude  that  there  is  a 
seasonal  variation  in  fallout.  The  present  discussion  supports 
this.  Note,  however,  that  there  are  anomalies  in  the  annual 
fallout  deposition  cycle,  in  that  fresh  debris  was  observed 
shortly  after  the  very  large  stratospheric  injections  of  1961 
and  1962. 
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3.  THE  PROBLEM  OF  NUCLEAR  CLOUD  RISE 
AND  DEBRIS  INJECTION  HEIGHT 


During  the  1961-1962  test  series,  some  38  thermonuclear 
devices  with  a  total  yield  of  approximately  3^0  Mt  were  fired, 
10%  in  tropical  and  90%  in  polar  latitudes  (Table  2-5).  For 
modeling  both  the  chemical  and  transport  effects  of  such  injec¬ 
tions,  one  needs  to  know,  or  at  least  to  hypothesize,  how  high 
the  clouds  went.  One  does  know  that  weapons  with  yields  greater 
than  several  megatons  deposit  most  of  their  debris  in  the 
stratosphere  and  that  the  effective  rise  height  increases  with 
increasing  yield. 

Figure  3-1  shows  several  models  of  cloud  rise  height  as 
a  function  of  yield.  The  wide  scatter  indicates  the  variabil¬ 
ity  and  uncertainty  of  the  results.  Thus  Peterson  (1970)  has 
constructed  a  model  of  cloud  rise  height  based  on  U.S.  data  in 
the  tropics  and  at  the  Nevada  Test  Site  (37°N)  for  low-yield 
tests.  Peterson  has  extended  his  model  to  high-latitude 
(Soviet,  75°N)  tests  based  on  considerations  of  atmospheric 
stability.  Note  that  the  largest  U.S.  test  was  a  15-Mt  sur¬ 
face  burst  (equivalent  airburst  yield  7.5  Mt),  so  beyond  this 
yield  there  are  no  experimental  data  to  support  his  model. 

Seitz  et  al.  (1968)  have  established  a  model  of  cloud  rise 
height  for  both  the  U.S.  and  USSR  1961-1962  test  series  as  an 

input  for  testing  the  transport  of  various  radioactive  tracers 

111  on 

and  thus  the  observed  profiles  for  C,  Sr,  and  other  tracers. 
See  Figs.  3-2  and  3-3  for  representative  profiles.  However, 
the  cloud  height  models  of  Seitz  et  al.  (1968)  for  the  Soviet 
bombs  shown  in  Fig.  3-1  are  not  the  result  of  direct  measure¬ 
ments.  Telegadas'  (1979b)  model  for  Chinese  (mid-latitude. 
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FIGURE  3-1.  Nuclear  cloud  rise  height  as  a  function  of  yield 
for  different  latitudes.  Peterson's  (1970)  equa¬ 
torial  model  comes  from  U.S.  tests  in  the  tropics 
and  (for  low  yields)  in  Nevada  (37°N).  Foley- 
Ruderman  (1973),  from  Chang  et  al .  (1979),  is  an 
ul trasimpl i f i ed  model  which  probably  gives  an 
upper  bound.  The  vertical  lines  labeled  P  or  T 
come  from  the  model  of  Seitz  et  al .  (1968);  the 
correction  for  height  above  the  tropopause  has 
been  made  by  Chang  et  al.  (1979)  and  provides  a 
plausible  lower  bound  for  cloud  rise.  The  verti¬ 
cal  line  labeled  "Mid-Latitude  (Telegadas)"  cor¬ 
responds  to  Chinese  3-Mt,  40°N  injections  (Tele¬ 
gadas  1977). 


3-2 


120,000 


100,000 


Km  20  P-* 1815  ' 
1  .1020 


•  2234' 

*  / 


80,000 


60,000 


I  5  r-.  n  ^5 


b  Ap-  AK*?PO?5r 


40,000 


20,000 


90°  60° 


30° 

NORTH 


30° 

SOUTH 


60°  90® 


FIGURE  3-2.  Atmospheric  burden  of  AHC,  December  1962--February 
1963,  mainly  from  aircraft  sampling.  (Source: 
Telegadas,  1971) 


approximately  3  Mt)  cloud  rise  is  based  mainly  on  aircraft  and 
balloon  observations  carried  out  mostly  along  the  west  coast  of 
North  and  South  America,  or  180°  in  longitude  from  the  point  of 
inj  ection  .. 

The  parameterization  of  Foley  &  Ruderman  (1973) — see,  e.g., 
Chang  et  al.  (1979) — has  also  been  used.  It  indicates  an  ex¬ 
tremely  high  cloud  rise  for  a  high-yield  bomb;  one  argument  for 
this  comes  from  the  observed  enhancement  of  the  lithium  air- 
glow  following  the  high-yield  Soviet  nuclear  explosions  of 
October  1961  (Sullivan  &  Hunten  1964).  If  the  enhanced  lithium 
airglow  extends  over  a  horizontal  dimension  of  the  order  of 
thousands  of  kilometers,  the  observed  airglow  levels  would  re¬ 
quire  several  grams  of  lithium.  Since  the  bomb  contained 
10-100  kg  lithium  (Sullivan  &  Hunten  1964),  this  observation 
is  hardly  convincing  evidence  that  a  significant  fraction  of 
the  cloud  rose  to  very  high  altitudes. 

In  view  of  this  uncertainty  in  injection  height,  I  shall 
follow  the  approach  of  Chang  et  al.  (1979),  who  use  two  dis¬ 
tinct  parameterizations — a  high  one  (Foley  Sc  Ruderman)  and  a 
low  one  (Seitz) — which  are  reduced  to  an  effective  mid-latitude 
injection,  as  Indicated  In  their  paper,  by  raising  the  cloud 
height  4  km  for  polar  tests  and  by  reducing  it  2  km  for  tropi¬ 
cal  tests.  These  two  parameterizations,  which  are  shown  in 
Fig.  3-1  (after  Chang  et  al.  1979),  represent  plausible  lower 
and  upper  bounds.  Using  these  two  models  for  debris  Injection 
height  and  dividing  the  1961-1962  Injections  into  five  distinct 
pieces,  each  with  an  appropriate  mean  yield  and  thus  effective 
injection  height,  gives  the  injection  groupings  listed  in 
Table  3-1.  The  major  stratospheric  injections  of  the  1950s 
are  grouped  comparably  in  Table  3-2. 

From  Table  3-1  one  can  infer  that  for  the  1961-1962  test 

14 

series  the  effective  injection  height  for  C,  which  is  a  tracer 
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TABLE  3-1.  GROUPING  OF  INJECTIONS  FOR  1961-1962  TEST  SERIES 


<0 


C£.  * 

t/l  t/1  •*  U 

(/)  XI  p  0) 
3E£Q(\J 

O  I  V£> 
X3  O  CD  03 
••  f—  3  r~ 

CM,—  l 


CM 

<£> 

03 


U 

O 


CO 

VO 

O 

OQ 

CM 

CM 

VO 

CM 

CM 

VO 

O 

O 

VO  CPv 
CM  r— 


CO  *>p  u 

(/>  to  r  oj 

ZD  -O  Q  CM 
EN  I  ^ 
o  cm  cnoi 

•*  -Cl  |  =3  r— 

r—  O  <C 
•  VO  00 
Q 


CO 

vO 

03 


u 

o 


ir> 

CM 


p 

r~  IO 

co 

oo 

o 


ir> 

vo 

to 

o 


vo 

CM 


vo  co 

CO  CM 


CM 

vo 

03 


CO  tO  p 
•  XD  u 
=>  E  O 

on| 


CM 

vo 

cr» 


3 

r? 


CO 


o 

o 


to 


CO 

VO 


vo  o% 

CM  r— 


P 

*  x 

or 

to  a>  *  co 

V/)  03  VI  to 

n  C  p 
E  «c 
o 

••  -o  to 
CO  CM  CMl 


vo 

03 


> 

o 


oo 

00 


VO 

CM 


CO 


no 

8 


vo 

# 

»— 


00 

co 


CM 


*  *  to 

OC  VI  Ol 
00  XI  |  r- 

to  E 

=>  O  CL  > 
X>  0>  O 
CO  z? 

••  co 

C  —  r-Kt 


CM  VO 
CM  »— 


C 

c 

o 

o 

p 

•r* 

u 

X 

V) 

p 

c 

VI 

o 

rtj 

* 

'r* 

P 

01 

TJ 

u. 

X 

•o 

c 

C  E 

01 

0) 

P 

» 

►-4  ^ 

XI 

•r- 

O 

XJ 

3 

>- 

p— • 

01  A 

or 

c 

01 

>  P 

i 

u 

C 

O  XI 

•r* 

mr~  -C 

>> 

N 

o 

•e— 

>» 

P  On 

0) 

P 

•r 

P  4) 

O 

•r- 

VI 

u  •*- 

C 

0>  01 

o 

0) 

VI 

<C  >- 

to 

V  X 

Lu 

VO 

V 

4- 

* 

p 

<V  = 
Er¬ 
rs  r— 

v)  ro 

vi  E 
ns  v) 

• 

0/  * 

01 

2:  0i 

u 

-c 

c 

•  <r  4J 

ro 

U 

VO  “T3 

*n 

vo  c 

4- 

03  <0 

c 

s _ _ 

CD 

0/ 

■r— 

1-  CD 

l/l 

01  4. 

>i 

P  <o 

r~ 

ro  = 

Q- 

oo 

E 

0> 

•r* 

P 

O 

-Q  P 

4-  P 

c 

Oi  P 

o 

-C  o 

p  p 

-o 

0) 

>1 

03  4. 

o  o 

(U 

p  p 

JC 

p 

CM  01 

</) 

VI 

XJ 

o 

C  Oi 

s*. 

ro  -C 

4. 

P 

01 

*  to 

P 

Q  -r- 

P 

O 

V)  "O 

XI 

CLr— 

c 

3  01 

3 

o  •*- 

o 

4^ 

4- 

O' 

01 

r— 

4.  ro 

o 

O  P 

3 

P  O 

XI 

P 

oi 

x» 

4. 

01  o 

p  p 

O 

o 

P 

3  XI 
CTr- 

>i 

01 

t— 

C  -r- 

c 

(V  >s 

o 

OJ 

JO  C 

X? 

O 

01 

to  T- 

p 

ro  V) 

o 

-C  to 

=3 

•r— 

cr 

X3  P 

oi 

0J  P 

u 

*r^  O 

ro 

>i 

O 

to 

C  •!“ 

01 

O  P 

4- 

•r-  (O 

3 

t0  4- 

cn 

lO 

•r» 

*»—  01 

P 

p  x: 

P  • 

01  *—  VI 

OJ  fO  P  P 


l-  P  «0 

o  -c 


2 


U 


3-6 


From  models  of  Chano  et  al .  (1979)— see  Fig.  3-1. 


TABLE  3-2.  GROUPING  OF  INJECTIONS  IN  THE  1950s 
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90 

of  total  yield,  was  some  1  to  2  km  greater  than  that  of  Sr, 
which  is  a  tracer  of  fission  yield,  mainly  because  the  Soviet 
58-Mt  bomb  of  30  October  1961  had  a  very  low  fractional  fission 
•yield . 

In  closing,  note  that  the  available  information  on  cloud 
rise  height  for  the  large  Soviet  bombs  of  1961-1962  is  based 
on  limited  and  indirect  data  and  should  be  used  with  extreme 
caution. 
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4.  EDDY  DIFFUSIVITY  PROFILES 


Figure  4-1  and  Table  4-1  show  some  eddy  diffusivity  or 
K-profiles  that  have  been  used  in  one-dimensional  (1-D)  model¬ 
ing  of  stratospheric  photochemistry  between  1974  and  the  present. 
The  following  points  should  be  noted: 

•  The  atmosphere  is  three-dimensional,  and  actual  trans¬ 
port  is  a  complex  combination  of  convective  and  dif¬ 
fusive  motions,  with  a  great  deal  of  space  and  time 
variability.  To  parameterize  this  complexity  in  terms 
of  down-gradient  diffusive  transport  with  simply  one 
parameter,  K  =  K(z),  which  depends  only  on  height,  is 
clearly  an  oversimplification.  The  1-D  parameterization 
has  been  used  extensively  in  describing  problems  which 
are  largely  determined  by  chemistry  but  where  the  dy¬ 
namics  must  be  inserted  in  the  simplest  feasible  way. 

It  is  relatively  satisfactory  for  processes  that  take 
place  sufficiently  far  above  the  local  tropopause,  per¬ 
haps  at  18  to  20  km,  for  details  in  the  region  of  the 
tropopause  to  be  unimportant.  In  the  application  to 
present-day  aircraft  Injections  which  often  take  place 
near  the  local  tropopause,  so  that  the  aircraft  some¬ 
times  fly  in  tropospheric  air  (high  moisture,  low  ozone) 
and  sometimes  in  stratospheric  air  (low  moisture,  high 
ozone),  the  concept  of  a  1-D  model  associated  with  a 
K-profile  must  be  used  with  extreme  caution. 

•  A  variety  of  different  K-profiles  have  been  suggested 
by  different  workers  (Fig.  4-1  and  Table  4-1).  They 
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0.263 
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FIGURE  4-1.  One-dimensional  eddy  diffusivity  profiles.  See 
Table  4-1  for  explanation. 
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all  show  similar  qualitative  characteristics — fast  trans¬ 
port  in  the  upper  stratosphere  and  in  the  troposphere, 
with  a  minimum  between  the  tropopause  and  altitudes  of 
18  to  25  km.  The  shape  and  magnitude  of  this  minimum  are 
critical  in  determining  how  fast  a  tracer  injected  in  the 
middle  or  upper  stratosphere  <>20  km)  is  swept  out  to  a 
sink  at  the  ground;  radioactive  tracer  data  provide  rela¬ 
tively  direct  information  on  this  minimum. 

•  Note  that  the  "slow"  profiles  such  as  WT  and  HN  have  a 
very  low  minimum  K-value  and  thus  provide  only  slow 
transport  out  of  the  stratosphere.  In  contrast,  the 
"fast"  profiles  DC  and  KD  have  a  relatively  large  mini¬ 
mum  and  so  produce  a  minimum  barrier  to  transport  out 
of  the  stratosphere. 

•  The  profile  "JA" — from  Jacobi  and  Andr£  (1963) — was 
devised  to  explain  the  behavior  of  fallout  products  in 
the  planetary  boundary  layer  (2  ^  1  to  2  km)  and  to  some 
extent  in  the  free  troposphere.  It  was  certainly  not 
Intended  to  be  used  in  the  stratosphere,  but  as  one  sees 
In  Fig.  4-1,  Its  behavior  below  20  km  is  not  Incon¬ 
sistent  with  some  other  profiles,  especially  DC  and  KD. 

•  The  rise  in  K-values  above  15  to  20  km  and  up  to  40  to 
50  km.  In  many  cases,  depends  on  the  N20  profiles  of 
Schmeltekopf  et  al.  (1977)  and  on  the  loss  mechanisms 
of  N20,  which  are  believed  to  be  approximately  80  per¬ 
cent  by  photolysis  with  solar  UV  below  230  nm,  and  20 
percent  by  reaction  with  0(^D),  about  half  of  which 
produces  NO.  The  older  estimates  of  K(z)  based  on 
methane  profiles  depend  on  the  assumed  loss  rate  by 
reaction  with  OH  and  0(^D)  in  a  chlorine-free  atmo¬ 
sphere;  the  loss  rate  varies  with  the  ambient  chlorine 
concentration,  among  other  factors.  There  is  also  the 
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fundamental  point  that  the  K-profile  depends  on  the  loss 
rate  of  the  tracer,  which  depends  on  the  model  of  photo¬ 
chemistry  for  a  reactive  tracer. 

•  The  most  recent  discussion  is  that  of  Danielsen  (1979), 
who  uses  three  very  different  tracers  (ozone,  nitrous 
oxide,  and  the  results  of  a  three-dimensional  computer 
simulation)  and  derives  a  profile  "KD",  which  in  his 
opinion  has  an  overall  uncertainty  of  a  factor  of  two. 

•  Except  in  the  critical  10-20  km  region,  few  of  the 
K-profiles  differ  by  more  than  a  factor  of  two. 

•  In  the  present  one-dimensional  analysis,  the  primary 
quantity  to  be  modeled  is  the  time-dependence  of  the 
stratospheric  burden  of  an  injectant.  For  this  I  use 
the  representative  description  of  Telegadas  (1971),  in 
which  the  stratosphere  is  taken  as  containing  20  per¬ 
cent  of  the  atmospheric  mass,  giving  a  global  mean  tro- 

popause  height  of  11  km.  Note  that  Ellsaesser  (1975), 

17 

p.  7-13,  suggests  a  stratospheric  mass  of  6.67  x  10 

-j  O 

kg  versus  Telegadas'  1.02  x  10  kg,  or  15  percent  of 
the  total  atmospheric  mass,  which  would  correspond  to 
a  global  mean  tropopause  height  of  approximately  13  km. 
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5.  THE  DIFFUSION  MODEL 

5.1.  MATHEMATICAL  FORMALISM 

The  atmosphere  is  characterized  by  a  number  density 
n  =  n(z)  which  depends  on  height  z;  under  isothermal  conditions 
(approximately  true  between  11  and  20  km),  one  has 

n(z )  =  n(zQ)  exp  [-(z-zq)/H],  (5.1) 

where  the  scale  height  H  =  kT/Mg  =  6.3  km  in  the  11-18  km 
region,  where  T  =  217°K.  Atmospheric  dynamics  is  parameter¬ 
ized  in  terms  of  the  K(z)  profile  discussed  in  Section  4, 
while  aerosols  are  characterized  by  a  sedimentation  speed 
u  =  u(z)  (see  Section  5.3). 

One  is  dealing  with  chemically  inert  tracers  characterized 
by  a  globally  averaged  mixing  ratio 

f /  t-v  _  burden  of  tracer  between  z  and  (z  +  Az)  at  time  t 

*  air  mass  between  z  and  (z  +  Az)  ’  ' 

f(z,t)  being  a  solution  of  the  following  diffusion  equation: 

(3/3z)  [n(z)K(z)3f/3z  +  n(z)u(z)f (z,t)]  =  n(z)3f/3t.  (5.3) 

The  differential  equation  (Eq.  5-3)  has  to  be  solved  for 
appropriate  initial  conditions  discussed  in  Section  5.2  and 
for  boundary  conditions  discussed  in  Section  5.4.  An  analytic 
solution  for  a  simplified  K-profile  was  given  earlier  (Bauer, 
Oliver,  Wasylkiwskyj ,  1978);  the  problem  is  here  solved  numeri¬ 
cally,  and  the  integration  scheme  is  outlined  in  Section  5.5 
and  discussed  in  detail  in  Bauer  (1980),  which  is  available 
on  request. 
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5.2.  INITIAL  CONDITIONS 


For  the  Chinese  explosions  it  is  assumed  that  the  mean  in¬ 
jection  height  is  approximately  18  km,  corresponding  to  the 
mid-latitude  injection  of  a  3-Mt  bomb.  Figure  3-4  shows  the 
vertical  Zr  activity  distribution  for  five  of  the  six  Chinese 
high-yield  explosions;  note  that  this  is  not  a  direct  observa¬ 
tion,  but  rather  has  been  reconstructed  largely  from  the  Air- 
stream  observations  flown  halfway  around  the  globe  from  the 
point  of  injection.  From  the  discussion  of  Telegadas  (1977), 
it  is  possible  to  assume  an  initial  distribution  for  the  mixing 
ratio  f(z,t)  of  Gaussian  form 

f(z,0)  =  fQ  exp  [-(z-z1)2/o^]  ,  (5.4) 

where  =  18  km  and  oq  =  2.15  km. 

The  model  of  Eq .  5.4  is  believed  to  be  adequate  for  a  mid¬ 
latitude  3-Mt  bomb,  which  corresponds  to  four  of  the  six  Chinese 
high-yield  tests.  The  other  two  Chinese  tests  had,  respectively, 
a  yield  of  "2  to  3  Mt"  (June  1973)  and  "4  Mt"  (November  1976) 

(see  Table  2-1),  and  to  a  first  approximation  one  may  again  use 
the  model  of  Eq.  5.4. 

For  the  Soviet  and  U.S.  atmospheric  nuclear  explosions  of 
1961  to  .1962,  for  which  the  yield  and  also  the  latitude  and 
season  of  injection  varied  significantly,  I  make  estimates  of 
injection  height  using  the  Foley-Ruderman  and  Seitz  models  of 
Section  3.  Note: 

•  For  a  given  injection  of  Table  3-1  or  3-2,  the  distri¬ 
bution  of  debris  mixing  ratio  is  described  by  a  Gaussian 
model  like  Eq.  5.4  with  mean  injection  height  zq,  the 
mean  of  the  lower  and  upper  bounds  of  total  activity 
from  Fig.  3-1,  and  with  oQ  =  2.15  km,  as  for  the  Chinese 
3-Mt  tests.  (The  choice  of  injection  heights  has  been 
discussed  in  Section  3;  I  use  both  "low"  (Seitz)  and 
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"high"  (Foley-Ruderman)  values.  It  would  be  possible 

to  use  different  values  of  o  .  such  as  one-fourth  of 

o 

the  spacing  of  lower  and  upper  bounds  of  Fig.  3-1, 

which  would  give  values  ranging  from  1.7  km  (Seitz, 

case  A  of  Table  3-1)  to  4.2  km  (Foley-Ruderman,  case  B 

of  Table  3-2).  Unlike  the  choice  of  injection  height, 

the  choice  of  a  is  not  critical,  since  we  are  asking 
°  gn  lli 

for  the  total  burden  of  Sr  and  C02  a  long  time 
after  the  injection,  and  thus  the  injection  height  and 
the  total  amount  of  material  injected  are  far  more  im¬ 
portant  than  the  details  of  the  initial  distribution 
function,  which  is  smeared  out  fairly  rapidly. 

•  For  the  U.S.  and  Soviet  injections  of  1961  and  1962, 
a  lumped  model,  replacing  the  28  distinct  explosions 
by  five  separate  injections,  is  developed  in  Section  3 
(see  Table  3-1  and  also  Table  3-2  for  stratospheric 
Injections  of  the  1950s).  In  effect,  a  series  of  in¬ 
jections  at  a  given  time,  location,  and  altitude  is 
combined  to  produce  a  single  averaged  Injection.  Refer¬ 
ring  to  Table  3-1,  group  A  corresponds  to  13  Soviet 
bombs  of  yield  In  the  1-  to  3-Mt  range  detonated  be¬ 
tween  1  September  1961  and  4  November  1961,  and  so  on. 
For  all  the  groupings  except  B,  the  ratio  of  fission 
yield/total  yield  Is  essentially  the  same,  1/3.  Group 
B,  consisting  of  two  very  large  Soviet  detonations 

(25  and  58  Mt )  had  a  very  small  ratio  of  fission  yield/ 

total  yield,  approximately  0.1,  and  as  a  result  of 

this  group — more  or  less  of  the  single  58-Mt  bomb — the 

1 4 

effective  injection  height  of  C  (which  is  proportional 

to  total  yield)  Is  approximately  1  to  2  km  higher  than 
90 

that  of  Sr,  which  is  proportional  to  fission  yield. 

14 

•  For  C  It  is  necessary  to  Include  the  earlier  nuclear 
explosions  to  simulate  the  behavior  after  1963.  The 
various  injections  from  1956  to  1958  are  simulated  by 
five  different  injections  listed  in  Table  3-2. 
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5.3.  SEDIMENTATION  OF  AEROSOLS 


90  95 

Sr  and  Zr  are  carried  on  stratospheric  aerosols,  sub¬ 
micron  particles  composed  mainly  of  ^SO^-HgO  mixtures  deposited 
on  various  condensation  nuclei.  These  aerosols  make  up  the 
"Junge  layer"  of  submicron  particles,  mainly  in  the  18-20  km 
lower  stratosphere  altitude  region.  The  particles  are  in  dy¬ 
namic  equilibrium;  photochemical  formation  of  from  sulfur 

compounds  diffusing  up  from  the  troposphere  (S02  and  COS,  mainly) 
and  reacting  with  OH  and  H02  radicals  in  the  stratosphere  is 
followed  by  a  variety  of  physical  processes  of  condensation, 
evaporation,  coagulation,  sedimentation,  and  diffusion.  The 
most  recent  and  comprehensive  model  of  the  Junge  layer  has  been 
given  by  Turco  et  al .  (1979)  and  Toon  et  al.  (1979);  see  also 
Hamill  et  al.  (1977),  who  discuss  the  physics  but  not  the  photo¬ 
chemical  production  mechanisms. 

The  average  of  a  surface-  and  volume-weighted  mean  of  the 
particle  size  distribution  computed  by  Turco  et  al.  (1979)  is 
0.3  ym,  but  Drevinsky  &  Peccl  (1965)  found  that  most  of  the 
radioactivity  is  carried  on  particles  of  0.02-0.15  ym  radius. 

Here  I  assume  a  rate  of  sedimentation  corresponding  to 
spherical  particles  of  radius  0.1  ym  (as  did  Telegadas  &  List, 
1969),  with  mean  density  of  approximately  1.8  g/cm^  (correspond¬ 
ing  to  saturated  H2S0i|/H20  droplets  at  these  altitudes).  I  do 
not  know  how  to  include  the  effects  of  all  the  physics  in  the 
recent  models  listed  at  the  beginning  of  this  section  in  a 
fallout  parameterization. 

A  natural  question  to  ask  is  how  significant  the  effect  of 
sedimentation  is  in  the  present  context.  Spherical  particles 
of  radius  0.1  ym  and  density  1.8  g/cmJ  have  a  sedimentation 
speed  u  =  0.0032  cm/sec  =  0.08  km/month  at  an  altitude  of  18  km 
(see  Hasten  1968,  whose  velocity  profile  differs  slightly  from 
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the  earlier  one  of  Junge  et  al. ,  1961).  For  injection  at  some 
height  z^  in  the  18  to  20  km  altitude  range,  an  aerosol  of  this 
size  is  transported  out  of  the  stratosphere  approximately  as 
fast  as  a  gas  injected  at  a  height  1-2  km  lower.  However,  for 
long-term  falloff,  the  effect  of  the  lower  boundary  condition, 
which  is  taken  up  in  the  next  section,  is  more  significant  than 
that  of  sedimentation. 

5.4.  BOUNDARY  CONDITIONS 

In  the  numerical  integration  of  the  diffusion  equation 
(Eq.  5.3)  one  has  injections  of  material  around  z  =  20  km, 
and  the  integration  extends  from  z  =  0  (ground)  to  z  =  50  km. 

The  choice  of  boundary  conditions  is  critical  for  the  inte¬ 
gration  and  differs  for  different  tracers.  At  the  upper  bound¬ 
ary,  the  proper  condition  to  choose  is  one  of  no  net  flux  across 
it.  The  flux  is  obtained  by  integrating  Eq.  5.3  with  respect 
to  z  for  the  flux  density: 

Net  flux  =  nK3f/9z  +  nuf.  (5.5) 

For  the  case  of  a  gas,  with  no  sedimentation  so  that  u  =  0,  the 
boundary  condition  of  no  net  flux  is 

9f/9z  =  0  at  upper  boundary.  (5.6) 

There  are  two  distinct  cases  for  the  lower  boundary  con¬ 
dition  . 

5.4.1.  Perfect  Sink  at  Bottom 

90 

This  condition  is  appropriate  for  fission  fragments  (  Sr 

OK 

and  Zr)  and  also  for  HTO  or  H^O  which  are  not  recycled  rap¬ 
idly  once  they  reach  the  ground,  and  it  is 

f  =  0  at  lower  boundary/surface.  (5.7) 

5.4.2.  Finite  Loss  Rate  at  Bottom 

14 

This  condition  is  appropriate  for  C02,  which  is  only 
absorbed  rather  slowly  by  the  ocean  (see,  e.g.,  Machta,  1973). 
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There  is  an  effective  loss  velocity  vg  characterized  as  follows: 

nK3f/3z  +  nuf  =  nvgf  at  lower  boundary.  (5.8) 

it 

The  value  of  v  for  C0~  is  obtained  as  follows.  Machta 
e  2 

(1973),  Pig.  2,  quotes  a  tropospheric  mass  of  C02  of  6  x  10  g, 
and  in  a  box  model,  a  transfer  rate  into  the  ocean  of  0.22 
year-1  (also  back  transfer  from  ocean  to  atmosphere  of  0.1 
year-1).  Now  the  concentration  of  li(C02  in  the  atmosphere 
after  the  nuclear  tests  was  far  greater  than  in  the  ocean,  and 
thus  back  transfer  should  be  neglected,  so  that  the  mass  trans- 
fer  rate  *  6  x  10  x  0.22  g/year  =  4.18  x  10  g/sec.  The 
density  of  CC>2  at  the  surface  is 

PCQ  =  1.225  x  10-3  x  3.3  x  10'^  g/cm3  =  4.04  x  10-7  g/cm3. 


Now  if  the  effective  surface  area  is  A,  then  the  deposition 
speed  vg  is  given  by  the  relation 


mass  transfer  rate 


(5.9) 


1 8  2 

Taking  A  =  total  surface  area  of  earth  =  5.11  x  10  cm  gives 

v  =  0.00202  cm/sec  *  0.053  km/month,  which  is  the  value  to  be 
e  14 

used  here  for  C02.  If  one  considers  that  most  of  the  exchange 
takes  place  with  the  ocean,  which  covers  two-thirds  of  the  sur¬ 
face  of  the  earth,  the  actual  exchange  velocity  with  the  ocean 
is  (3/2)  x  0.053  km/month  =  0.08  km/month.  However,  this  is 
not  the  appropriate  velocity  to  be  used  here,  since  we  are 
dealing  with  a  global  mean  velocity. 


5.5.  THE  NUMERICAL  MODEL 


The  left-hand  side  of  Eq.  5.3  contains  two  terms: 


(3/3z ) 

(nK3f/3z ) 

(a) 

(3/3z) 

(nuf). 

(b) 

(5.10) 
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Because  here  is  more  important  than  T^,  the  problem  is  the 
integration  of  a  parabolic  partial  differential  equation,  for 
which  the  ordinary  explicit  numerical  integration  scheme  be¬ 
comes  unstable  unless 

2  KAt/(Az)2  «  1.  (5.11) 

It  is  appropriate  to  use  spatial  steps  of  Az  =  1  km  (note  also 

that  the  K-profiles  that  are  normally  used  have  discontinuities, 

which  may  cause  further  difficulties),  and  because  the  actual 

1}  2  2 

K-values  are  of  the  order  of  10  cm  /sec  (2.63  km  /month)  in 
the  lower  stratosphere,  or  103  cm  /sec  (26.3  km  /month)  in  the 
troposphere  and  upper  stratosphere,  the  condition  of  Eq.  5.11 
requires  a  time  step  as  small  as  At  =  0.01  month  for  stability 
if  one  uses  an  explicit  scheme.  With  the  implicit  Crank- 
Nicholson  scheme,  which  is  customarily  used  for  parabolic 
equations  of  this  type  (see,  e.g.,  McCracken  &  Dorn,  1964, 
pp.  365-404;  Richtmyer  &  Morton,  1967,  pp .  185-201),  a  time 
step  At  =  1  month  gives  stable  and  accurate  results,  as  checked 
against  the  analytic  two-step  solution  of  Wasylkiwskyj  (see 
Bauer,  Oliver,  &  Wasylkiwskyj ,  1978,  and  Bauer,  1980,  item  3). 
With  a  time  step  At  =  1  month,  running  for  times  t  =  20  months, 
the  program  takes  several  seconds  to  run  on  IDA's  CDC-6400. 

Times  as  long  as  t  =  15  years  are  needed  for  simulation  of  the 
1961-1962  test  series,  and  thus  it  is  evident  that  an  explicit 
scheme,  with  a  time  step  At  =  0.01  month,  would  be  extremely 
time  consuming  and  subject  to  large  roundoff  errors  in  addition 
to  possible  errors  associated  with  discontinuities  in  the  K- 
profile.  The  implicit  Crank-Nicholson  scheme  is  stable  and 
accurate.  A  complete  documentation  of  the  computer  program 
is  given  in  Bauer  (1980),  which  is  available  on  request. 

5.6.  IMPULSIVE  VERSUS  STEADY-STATE  INJECTIONS:  A  COMPARISON 
OF  ATMOSPHERIC  RESIDENCE  TIMES 

One  may  ask  for  the  relation  between  the  present  model  for 
an  impulsive  injection  and  the  case  of  a  steady-state  strato- 
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spheric  injection,  such  as  that  due  to  aircraft.  For  a  steady- 
state  injection,  one  may  define  an  atmospheric  residence  time 

T  =  flux/burden  (5.12) 

s  s 


(see,  e.g.,  Bauer  &  Gardner,  1977).  For  the  present  case  one 
may  define  a  residence  time  T^  in  terms  of  the  time-dependent 
stratospheric  burden  B  ,  , 


Bstrat(tl  + 


T,  )  = 
imp 


(1/e>  Bstrat(tl> 


(5.13) 


i.e.,  as  that  time  in  which  the  stratospheric  burden  falls  to 
1/e  of  its  value  at  some  initial  time  t,. 

Some  numerical  examples  for  an  18-km  injection  and  a  DC 

profile  are  given  in  Table  5-1,  which  shows  that  for  this  case 

there  is  a  significant  variation  in  Tlmp  as  t^  changes.  (Note 

that  for  a  "slow"  profile  like  WT  much  less  variation  in  T, 

imp 

would  be  detectable  over  the  50-month  run  from  which  the  esti¬ 
mates  of  Table  5-1  were  made.) 


TABLE  5-1.  STRATOSPHERIC  RESIDENCE  TIMES  FOR  STEADY 
STATE  AND  IMPULSIVE  INJECTIONS 


(K-proflle:  DC,  18-km  injection  of  a  gas,  no 
recycling  at  ground) 

Steady  State,  Eq.  5.12: 

Tss  =  18.9  months  (Bauer  and  Gardner,  1  977  ,  p.  D - 1 9 ) 
Impulsive  Case,  Eq.  5.13: 

Initial  Time  Residence  Time, 

t, ,  months  T .  ,  months 

_J _  imp  _ 


0 

12.1 

5 

13.6 

10 

18.9 

15 

23.8 

20 

27.2 
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In  conclusion,  it  should  be  noted  that  while  one  K-profile 
may  be  applied  to  a  steady  or  impulsive  injection  at  arbitrary 
heights  and  be  used  to  describe  the  parameterized  dynamics,  yet 
an  atmospheric  residence  time  can  be  defined  in  different  ways 
which  may  give  different  numerical  values  even  if  the  underly¬ 
ing  dynamics,  i.e.,  the  K-profile,  are  the  same.* 


it 

For  a  detailed  discussion  of  various  definitions  of  residence 
times  under  steady  and  unsteady  state  conditions,  see,  e.g., 
Schwartz  (1979). 
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6.  COMPARISON  OF  THE  DATA  WITH  CALCULATIONS 


6.1.  INTRODUCTION 

Section  6.2  compares  the  decline  in  stratospheric  burdens 
of  -^Zr  from  several  Chinese  thermonuclear  explosions  with 
model  predictions.  Section  6.3  treats  the  decline  in  strato¬ 
spheric  burden  of  ^°Sr  following  the  1961-1962  injections, 

lii 

while  Section  6.H  considers  the  time  variation  in  excess  C 0o 

ill  ^ 

through  1969;  because  of  the  slow  loss  rate  of  C02  to  the 
ground  (mainly  into  the  ocean),  it  is  necessary  to  consider 
injections  from  the  1950s  (Table  3-2)  as  well  as  the  larger 
injections  from  1961-1962  (Table  3-1). 

These  analyses  of  the  decay  in  stratospheric  burden  with 
time  require  that  the  initial  injection  height  of  the  tracers 
should  be  known.  Especially  for  the  Soviet  bombs  (which  con¬ 
tributed  90%  of  total  yield  during  the  1961-1962  time  period) 
the  injection  heights  are  not  well  known  (see  Section  3  for  a 
discussion  of  the  problem),  and  thus  in  Section  6.5  I  adopt 
an  alternative  approach  which  was  first  used  by  Johnston, 

Kattenhorn  &  Whitten  (1976),  who  analyzed  the  time  variation 

on  14 

of  the  altitude  profiles  of  Sr  and  C02.  My  conclusions 
are  rather  different  from  those  of  Johnston  et  al.,  and  this 
is  traced  back  to  the  boundary  conditions  used. 

6.2.  CHINESE  BOMBS:  95Zr  DATA,  1967-1978 

The  discussion  of  Section  2.3,  and  in  particular  Pig.  2-6, 

95 

shows  that  when  the  decay-corrected  stratospheric  burdens  of 
are  normalized  and  compared  on  a  seasonally  adjusted  basis, 
they  show  slow  fallout  until  the  first  winter  season  and  fairly 
rapid  fallout  thereafter. 
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Figure  6-1  compares  calculations  with  the  observed  and 
or 

suitably  normalized  Zr  data.  The  calculations  all  assume  a 
mean  injection  at  18  km  (see  Eq.  5*4  for  the  initial  condition  as¬ 
sumed)  and  the  following  boundary  conditions  for  the  mixing 
ratio  f(z,t)  (Eq.  5.2): 

Perfect  sink  at  the  bottom,  i.e.,  f(z  =  0,t)  =  0  (6.1) 

x’lux  through  the  top  (Section  5.4),  (6.2) 

which  are  appropriate  for  this  precipitation  scavengeable 
tracer. 

I  assume  sedimentation  (see  Section  5*3  for  the  model 
used),  and  injections  on  1  January  and  1  July  for  fast  trans¬ 
port  (Kp^,)  and  slow  transport  (K^),  respectively. 

Overall,  fast  transport  for  a  January  injection,  or  after 
a  delay  until  the  first  winter  following  a  summer/fall  injec¬ 
tion,  provides  an  adequate  representation  of  the  data. 

G . 3  .  90Sr  DATA,  1  960-1  967 

If  one  looks  at  Fig.  2-7,  one  sees  that  in  comparison  with 
1  h  on 

C  the  stratospheric  burden  of  Sr  decays  rather  rapidly 

with  time,  and  so  in  analyzing  the  data  it  is  not  necessary  to 
correct  for  injections  during  the  1950s.  The  injection  of  the 
1961  to  1962  test  series  is  modeled  in  terms  of  the  five  dis¬ 
tinct  injections  shown  in  Table  3-1.  Figure  6-2  compares 

90 

the  stratospheric  Sr  burden  from  Fig.  2-5  with  model  calcu¬ 
lations  using  the  DC  and  WT  profiles  of  Section  4  (with  the 
effects  of  sedimentation  Included),  and  with  the  "Foley-Ruderman" 
(hie-h)  and  "Seitz"  (low)  cloud  rise  estimates  as  used  by  Chang 
et  al.  (1979).  The  absolute  magnitude  of  the  computed  burden 
is  fitted  to  the  observed  activity  in  January  1963. 

One  sees  that  the  fast  DC  transport  with  the  low  (S)  cloud 

90 

rise  height  fits  the  decay  of  Sr  after  January  1963  far  better 
than  any  other  combination.  Fast  transport  is  consistent  with 
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95 

the  Zr  results,  since  the  Soviet  injections  took  place  in 
late  fall,  so  that  the  fast  winter  transport  took  over  shortly 
after  the  injection.  Regarding  the  cloud  rise,  evidently  the 
lower  the  injection  height,  the  faster  the  removal  for  a  given 
K-profile,  and  certainly  the  Airstream  data  offers  no  evidence 
for  the  PR  (high)  cloud  rise. 

After  1967,  the  effects  of  Chinese  and  French  nuclear  ex- 

90 

plosions  raised  the  Sr  level  to  a  detectable  extent,  so  the 
present  comparison  terminates  there. 

6.4.  14C  DATA,  1958-1969 

90 

Reference  to  Fig.  2-7  shows  that  in  comparison  with  Sr, 

which  falls  out  of  the  stratosphere  rather  rapidly  and  is  rained 

out  of  the  troposphere  so  that  its  tropospheric  burden  is  very 
14 

small,  C,  which  is  found  in  the  atmosphere  as  C0~,  stays  much 

14  d 

longer.  While  the  higher  injection  height  of  C  and  the  sedi- 
90 

mentation  of  Sr  do  have  an  effect,  the  principal  cause  of 

14 

this  has  to  do  with  the  lower  boundary  condition:  C02  is 

only  lost  very  slowly  at  the  surface — mean  disappearance  time 

4.5  years  (see  Maehta  1973) — and  is  present  in  the  troposphere 

for  a  number  of  years  after  stratospheric  injection  (Fig.  2-9). 

14 

Thus  the  lower  boundary  condition  for  CC>2  is  quite  different 
from  that  for  any  of  the  other  tracers  considered  here.  See 
Section  5.4.2  for  a  discussion. 

14 

Because  of  the  long  atmospheric  residence  time  of  C0o, 

14  d 

if  one  wishes  to  study  the  effects  of  C  injections  in  1961- 
1962,  it  is  necessary  to  consider  the  injections  of  the  1950s 
(Table  3-2)  as  well  as  the  1961-1962  injections,  which  are 
parameterized  as  shown  in  Table  3-1.  The  earlier  injections 
provide  a  background  for  the  1961-1962  tests,  and  thus  I  com¬ 
pare  results  between  1958  and  1969,  after  which  the  effects  of 
French  and  Chinese  explosions  obscure  the  record. 
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The  data  of  Fig.  2-6  over  the  period  1958-1969  are  thus 
compared  with  the  results  of  the  present  simulation,  using  the 
high  (FR)  and  low  (S)  injection  heights  of  Section  3  and  the 
"fast"  (DC)-  and  "slow"  (WT)  K-profiles  of  Section  4 .  Figure 
6-3a  shows  calculations  for  high  injection  height  ("Foley- 
Ruderman"  of  Section  3»  or  FR) ,  and  Fig.  6-3b  shows  the  results 
for  low  injection  height  ("Seitz"  of  Section  3,  or  S).  One  sees 
that : 

•  The  high-altitude  injection  gives  too  slow  a  decrease 

i  h 

in  the  stratospheric  burden  of  excess  C  and  also  too 
late  and  too  small  a  maximum  in  the  tropospheric  burden. 

•  The  low  injection  height  with  a  "fast"  transport  co¬ 
efficient  reproduces  both  the  observed  decrease  in 

stratospheric  burden  and  also  the  peaking  in  tropo- 
14 

spheric  C. 

6.5.  ANALYSIS  OF  ALTITUDE  PROFILES  OF  14C:  DIFFERENCES 
BETWEEN  THE  ANALYSIS  OF  JOHNSTON,  KATTENHORN,  AND 
WHITTEN  (1976)  AND  THE  PRESENT  WORK 

Johnston,  Kattenhorn  &  Whitten  (1976)  took  the  altitude- 
latitude  network  of  excess  data  from  1963  to  1969  (Telegadas 
1971)  and  constructed  a  Northern  Hemisphere  mean  altitude  pro¬ 
file  of  both  number  density  and  mixing  ratio  at  various  times. 
They  found  that  this  Northern  Hemisphere  mean  (corrected  for 
leakage  to  the  Southern  Hemisphere)  agrees  well  with  the  ob¬ 
served  profile  from  30°N,  so  that  in  this  instance  two  repre¬ 
sentations  used  for  a  1-D  parameterization  agree.  They  used 
a  numerical  diffusion  equation  to  test  how  well  different  K- 

profiles  were  able  to  reproduce  the  changing  altitude  profile 
14 

of  C,  mainly  on  a  year-by-year  basis,  but  also  over  longer 
times;  and  they  found  that  the  Hunten  (HN)  profile,  a  fairly 
slow  K-profile,  reproduces  most  of  the  data  better  than  any  of 
the  other  K-profiles  they  considered. 

Their  boundary  conditions  are  described  in  their  paper  as 
follows:  "The  lower  boundary  condition  was  that  observed  at 
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excess  C  in  stratosphere  and  troposphere 
ection  height. 


1  km,  which  remained  constant  for  several  years  near  3  x  10“^ 
mixing  ratio.  The  upper  boundary  condition  was  that  the  con¬ 
centration  at  51  km  was  one  half  that  at  50  km." 

li| 

These  boundary  conditions  are  not  appropriate  for  C  or 
1 24 

C02,  because  they  are  equivalent  to  very  strong  sinks  at  both 
the  upper  and  lower  boundaries. 

Consider  first  the  lower  boundary.  In  the  present  nota¬ 
tion,  the  mixing  ratio  at  the  bottom  step  is  f,  =  c  =  constant 

—1 6  ^ 

(3  x  10  D).  This  means  that  there  is  a  very  strong  sink  at 

z  =  1  km  because  any  material  that  comes  to  the  bottom  step, 

tending  to  change  f  from  the  value  c,  is  immediately  absorbed. 

Formally,  the  mixing  ratio  f(z,t)  is  a  solution  of  a  linear 

(diffusion)  equation,  Eq.  5.3,  which  one  may  write 

L(f)  =  0.  (6.3) 

Now,  let 

f(z,t)  =  c  +  g(z,t).  (6.4) 

Then,  since  c  is  a  constant,  L(c)  =  0,  and  thus  g  (z,t)  is  a 
solution  of  the  diffusion  equation,  one  has  L(g)  =  0.  Thus  an 
equivalent  way  of  describing  what  is  done  is  to  solve  the  dif¬ 
fusion  equation  L(g)  =  0  with  the  lower  boundary  condition 

gx  =  fx  -  c  =  0  (6.5) 

and  then  to  add  the  constant  c  to  the  profile  at  all  altitudes. 

The  boundary  condition  of  Eq.  6.5  is  appropriate  for  a  tracer 
90  95 

such  as  Sr  or  Zr  which  is  precipitation  scavenged  and  lost 

14 

on  the  ground,  but  not  for  C02,  for  which  the  discussion  of 
Section  5.4.2  shows  that  a  much  weaker  sink  exists  at  the 
bottom:  explicitly,  Eq.  5.8  shows  that  the  lower  boundary  con¬ 

dition  is 

9f/3z  =  vef/Kb  ,  (6.6) 

p 

where  vg  *  0.053  km/month,  and  ~  25  km  /month. 
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The  upper  boundary  condition 

f  =  —  f 
51  2  1 50 

defines  a  net  sink  at  the  top,  because 

grad  f  =  3f/3z  =  (f^-f^J/CAz)  =  -  f5Q/2(Az)  -  -  f/K,.,  (6.8) 

where  now  v  =  K./2(Az)  ~  10  km/month.  It  can  readily  be  shown 
(see,  e.g.,  item  3  of  Bauer  1980)  that  this  boundary  condition 
corresponds  to  a  very  strong  sink,  approximately  equivalent  to 
a  "perfect  sink  "  f50  -  (The  transformation  f  =  g  +  c  changes 
the  numerics  a  little,  but  not  the  essence  of  the  argument.) 

As  a  consequence  of  these  boundary  conditions,  Johnston 

et  al.  (1976)  have  a  source  of  at  approximately  20  km  and 

study  the  loss  of  material  to  sinks  at  the  ground  (1  km)  and 

also  at  the  stratopause  (50  km).  Because  these  sinks  are  un- 

14 

physically  strong  for  CC^,  they  must  have  unphysically  slow 
transport  between  the  source  and  the  sinks  to  reproduce  the  ob¬ 
served  concentration  profile. 

Johnston  et  al.  (1976)  reinitialize  their  concentration 
profile  at  annual  intervals,  which  reduces  the  effects  of  the 
boundary  conditions  and  of  leakage  to  the. Southern  Hemisphere, 
but  which  can  lead  to  Incorrect  conclusions  if  there  is  a 
change  in  effective  K-profile  with  time  such  as  was  found  by 
Mahlman  (1975)  in  a  3-D  numerical  simulation. 

I  have  done  calculations  similar  to  those  of  Johnston 
et  al.,  using  their  initial  profile  of  January  1963  (see  their 
Table  5)  in  the  form 

fic(z)  =  fcoexp[-(z-zc)2/0c2]»  (6-9) 

where  f  =  91.5xl0”1^,  z  *  21.4  km,  a  ■  3.65  km,  and  all 

CO  c  c 

permutations  of  the  following  boundary  conditions  are  used: 
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Upper  boundary: 

(a) 

no  flux,  8f/8z  =  0 

(b) 

perfect  sink,  f  *  0 

Lower  boundary: 

(a) 

proper  boundary  condition  for 
i.e.,  Eq.  6.6 

(b) 

perfect  sink,  f  =  0. 

Note  that  these  conditions  are  not  identical  with  those  of 
Johnston  et  al.,  but  they  are  very  similar. 

Some  representative  mixing  ratio  profiles  are  shown  in 
Pig.  6-4  for  January  1965,  i.e.,  two  years  after  the  initial 
value  of  Eq.  6.9.  The  results  are  shown  for  the  (relatively 
slow)  HN  profile,  and  one  sees  that  the  effects  of  both  boundaries 
have  traveled  to  the  maximum  in  the  mixing  ratio.  The  faster 
DC  profile  gives  a  somewhat  lower  maximum  mixing  ratio,  of  the 
order  of  20-25xl0-1^ ,  which  is  in  better  agreement  with  experi¬ 
ment.  However,  at  this  time  the  DC  mixing  ratio  does  not  de¬ 
crease  above  30  km,  which  suggests  that  perhaps  the  high-alti¬ 
tude  transport  given  by  this  model  is  too  fast,  and  that  the 
high-altitude  value  of  the  HN  profile  is  more  appropriate. 

The  faster  DC  profile  also  leads  to  more  rapid  transport  of 
the  boundary  conditions  toward  the  maximum. 

Figure  6-5  compares  the  observed  global  and  Northern 
Hemisphere  stratospheric  burdens  of  ^C  with  these  computations 
(i.e.,  the  initial  condition  of  Eq.  6.9)  for  January  1963  with 

14 

DC  and  HN  K-profiles,  and  for  the  proper  CO 2  boundary  con¬ 

ditions  as  well  as  for  the  case  of  a  perfect  sink  at  both  top 
(50  km)  and  bottom.  One  sees  that,  as  would  be  expected,  the 
effect  of  rapid  transport  with  weak  sinks  tends  to  compensate 
for  the  case  of  slow  transport  with  strong  sinks.  Note  that 
the  data  provide  a  better  fit  to  the  global  than  to  the  Nor¬ 
thern  Hemispheric  burden,  the  difference  being  transport  to 
the  Southern  Hemisphere. 
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FIGURE  6-5.  Decrease  in  burden  of  i_rC02  following  the  initial 
value  of  Eq.  6.9.  Note  effects  of  varying  K- 
profiles  and  boundary  conditions  in  the  modeling, 
and  using  global  or  Northern  Hemisphere  data  (see 
discussion  In  Section  6.5).  "Sinks"  implies 
perfect  sinks,  f  ■  0,  at  both  upper  and  lower 
boundaries . 


In  conclusion,  the  present  calculation  differs  from  the 
work  of  Johnston  et  al.  (1976)  in  considering  global  strato¬ 
spheric  and  tropospheric  burdens  rather  than  hemispheric  al¬ 
titude  profiles,  in  emphasizing  longer  times  after  injection, 

and  in  a  better  treatment  of  the  boundary  conditions,  so  that 
on  14 

the  behavior  of  Sr  can  be  explained  as  well  as  that  of  C02. 
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7.  CONCLUSIONS 


The  following  conclusions  may  be  drawn  from  the  analysis 
of  this  paper: 

95 

1.  The  Zr  stratospheric  data  from  six  Chinese  thermo¬ 
nuclear  explosions  demonstrate  a  strong  seasonal  vari¬ 
ation  in  the  decrease  in  stratospheric  burden,  with 
little  stratosphere-to-troposphere  transfer  for  a 
summer  or  fall  injection  until  the  next  winter  sea¬ 
son.  For  a  winter  injection,  rapid  stratosphere-to- 
troposphere  transfer  starts  immediately.  (There 

have  been  no  Chinese  explosions  in  springtime.) 

2.  As  far  as  decay  in  stratospheric  burden  is  concerned, 
these  data  can  be  modeled  in  terms  of  one-dimensional 
(1-D)  transport  with  a  ’'fast"  profile  (such  as  DC 

or  KD  of  Section  Fig.  4-1)  for  a  winter  injection. 
For  a  summer  or  fall  injection  it  is  again  best  to 
use  a  "fast"  K-profile,  but  with  the  time  of  injec¬ 
tion  delayed  until  the  next  winter. 

90  14 

3.  The  Sr  and  excess  C  stratospheric  burdens  due  to 

the  USSR  and  U.S.  nuclear  tests  of  1961-1962  have 

been  analyzed  with  an  extension  of  the  same  technique 

to  take  into  account  the  variation  in  yield  and  lati- 

90 

tude  of  the  different  bombs.  The  Sr  falls  out  much 

ill 

more  rapidly  than  does  the  C,  principally  because 

90 

of  the  difference  in  lower  boundary  conditions.  Sr 

in 

is  rained  out  rapidly  and  completely,  while  C  (which 
is  carried  mainly  as  ^COg)  is  only  absorbed  slowly 
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(e-folding  time  of  4.5  years)  by  the  ocean  (and  bio¬ 
sphere).  Additional  contributions  to  the  difference 

90 

come  from  the  sedimentation  of  Sr  and  from  its  lower 
effective  injection  height. 


an 

The  falloff  in  stratospheric  burden  of  Sr,  and  of 

1 4 

stratospheric  and  tropospheric  burdens  of  C,  can 
be  fitted  well  by  using  a  "low"  injection  height 
("Seitz"  model  of  Section  3)  and  a  "fast"  K-profile 
(such  as  DC  or  KD  of  Section  4).  Any  other  combina¬ 
tion  of  injection  height  and  K-profile  gives  signifi¬ 
cantly  slower  decay  in  burdens  of  these  tracers  than 
the  observations  show. 


5.  Reference  to  Pig.  2-7  shows  that  the  inferred  strato¬ 
spheric  burden  of  HTO  falls  off  somewhat  more  slowly 

90 

than  that  of  Sr  and  (especially  at  later  times) 

14 

significantly  faster  than  that  of  CO^.  This  can 
be  explained  as  follows: 

a.  T  is  a  tracer  of  fusion  rather  than  fission,  so 

its  effective  injection  height  is  somewhat  lar- 

14 

ger  than  that  of  C  and  definitely  larger  than 
that  of  ^°Sr. 

b.  HTO  is  a  gas,  so  it  is  not  subject  to  sedimenta- 

90 

tion  as  Sr  is. 

c.  HTO  is  precipitation  scavenged  in  the  troposphere 
and  is  lost  on  the  ground  after  a  large  injection, 
so  its  lower  boundary  condition  resembles  that  of 
^Sr  rather  than  that  of  ^COg. 

6.  It  is  not  necessary  to  postulate  extra-rapid  sedi- 

90 

mentation  of  7  Sr  or  a  large  unobserved  high-altitude 
14 

reservoir  of  C  to  explain  the  difference  in  decay 
of  these  tracers. 
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7.  The  interpretation  of  Johnston  et  al.  (1976)  that  the 
profile  data  require  relatively  slow  (HN)  trans¬ 
port  is  shown  to  be  related  both  to  their  selection 
of  relatively  short  sampling  intervals  and  to  their 
incorrect  choice  of  boundary  conditions. 
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